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Abstract. Local-box hydrodynamical model atmospheres provide statistical information about the spatial dependence, as well 
as temporal evolution, of a star's emergent radiation field. Here, we consider late-type stellar atmospheres for which temporal 
changes of the radiative output are primarily related to convective (granular) surface flows. We derived relations for evaluating 
the granulation-induced, disk-integrated thus observable fluctuations of the stellar brightness and location of the photocenter 
from radiation intensities available from a local model. Apart from their application in the context of hydrodynamical stellar 
atmospheres, these formulae provide some broader insight into the nature of the fluctuations under consideration. Brightness 
fluctuations scale inversely proportional to the square root of the number of convective cells (the statistically independently 
radiating surface elements) present on the stellar surface and increase with more pronounced limb-darkening. Fluctuations of 
the stellar photocentric position do not depend on the number of cells and are largely insensitive to the degree of limb-darkening. 
They amount to a small fraction of the typical cell size, and can become a limiting factor for high-precision astrometry in the 
case of extreme giants. The temporal brightness and positional fluctuations are statistically uncorrelated but closely related in 
magnitude. 
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1 . Introduction variability would aid the planning and interpretation of upcom- 

. £^ ' ing observations. 

^ ! The presence of granular flows on the surfaces of late-type Rased Qn detailed radiation . hydrodynamics simulations of 

^ ; stars implies ultimate limits to their photometric stability, the convect i ve surfac e flows - similar to the ones described in 

stability of the position of their photometric centroid (here- | Wedemever et al ] £ 20 04) - we want to predict the properties of 

after "photocenter"), and the stability of their spectroscopi- granu , ation _ re i ated micr o-variability in late-type stars. To our 

cally determined radial velocities. The searches for extrasolar knowledge the work of IrramDedachttli] &) is the only 

planets and stellar p-mode oscillations have been pushing to- example of such a meoretical effort t0 date T rampedach et al. 

wards higher observational sensitivity, now allowing the de- smdied the brightness and radial velocity variability in the Sun, 

tection of the low-level (hence "micro-") variability induced ff Cen A> and Procyon A Their hydrodynamical models> as 

by the stochastical ly changing granular flow pattern (e.g., weU afj ^ Qnes wfi refgr ^ ^ so . called - local _ box » models . 

| Kieldsen et al, | 1999|) . Stellar micro-variability becomes partic- They simulate convection in a smalL representative volume lo- 

ularly important for present and future space missions such as cated wkhin the stellar surface layers In order tQ derive obsery . 

Most, Corot, Kepler, Sim, and Gaia (|Hatzes|2002HGreen et alj ^ disk . integrated quantities, properties of the local model 



O 

in 
o 

Oh 
6 



|200J |AigrainetalJ | 2004t | Matthews et al. || 2004 | h In this con- haye tQ be extrapolated to the full stellar disk 

text, granulation-related micro-variability is usually considered _ . ,. • • , , 

. _ i iii • ,, , ,• This paper outlines the statistical procedure that we employ 

as a noise source. It should be appreciated, however, that this „ , T .. ,, , , 

. „ , . . • , , • • r for that purpose. It is intended as a precursor to later work pre- 

noise also carries information about the statistics of con- . ... , ,,„.,; • • , , 

. • , , ■ • senting simulation results. While the statistical developments 

vective flows. Whether considered as noise or signal: a bet- , . , _ , 

, • , , r , ■ r t • are mostly straightforward, we nevertheless feel that a some- 

ter theoretical characterization ot granulation- related micro- . , 

what detailed discussion is necessary to judge later results and 

- perhaps even more important - their inherent uncertainties. 

Send offprint requests to: Hans-Giinter Ludwig In the following we consider brightness fluctuations and flue- 
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tuations in the position of the stellar photocenter. We do not 
consider the important case of fluctuations in spectroscopic ra- 
dial velocity, basically since related spectral line formation is- 
sues complicate the problem. Despite its mostly technical na- 
ture, this paper contains some results within a broader scope 
so that the non-technically oriented reader might also find it of 
interest. 

The paper is organized in the following way: in Sect. [2] we 
start by stating the basic statistical assumptions made when ex- 
trapolating from the local simulations to the full stellar disk, 
and describe in Sect. [3] what kind of simulation data did en- 
ter. In Sect. @]we consider statistical properties of an individ- 
ual simulation patch, and in Sect. [5] we continue by deriving 
the related formulae for the expectation value, variance, and 
power spectra of temporal brightness and positional variations 
of disk-integrated quantities. We show that there is a general 
close connection between brightness and positional variations, 
and that - a somewhat counter-intuitive result - photocentric 
variations do not depend on stellar radius but on the atmo- 
spheric parameters alone. Along the way, we introduce an ap- 
proximation related to the well-known Eddington approxima- 
tion to mitigate some numerical shortcomings. In Sect. |6] we 
add some cautionary comments about short-cuts when evaluat- 
ing supposedly disk-integrated quantities from the flux data of 
a simulation box. In Sect. we discuss an illustrative analyti- 
cal model of brightness and positional variations, and make the 
connection between simulation patches and convective cells. In 
Sect.|S]we give our conclusions. 

2. Basic statistical assumptions 

Our goal was to make predictions of a star's brightness fluctu- 
ations and of the displacements of its photocenter introduced 
by granulation, resorting to results of radiation-hydrodynamics 
simulations. These simulations provide a statistically represen- 
tative, small, rectangular patch (or "tile") of the emergent radi- 
ation field and its temporal evolution on the stellar surface. We 
imagined the visible stellar hemisphere as tiled by a possibly 
large number of simulation patches. Luckily, we did not have 
to construct this tiling, which would have been difficult in prac- 
tice due to the periodic lateral boundary conditions employed 
in the simulations, among other problems. To derive what a dis- 
tant observer would record if observing the unresolved stellar 
disk, we made the following statistical and physical assump- 
tions: 

- Each patch radiates statistically independent of all other 
patches. 

- All patches on the stellar surface share the same statistical 
properties. 

- The space-wise statistics of the patches on the stellar sur- 
face can be replaced by the time-wise statistics of a simu- 
lated patch. 

In a real star the first assumption is not exactly fulfilled. In 
particular, acoustic modes have correlation lengths which in- 
troduce a coupling over scales larger than are typically encom- 
passed by a simulation box. The coupling between modes and 
convection is too weak to introduce a mode-related long-range 



correlation into the convective part of the flow, so that the first 
assumption should hold for granulation itself. The modal part 
could be treated separately as a perfectly correlated signal. In 
supergiants one might encounter limits to a statistical treatment 
of convection-related variability, since only very few granules 
cover the stellar surface (e.g., Frevtas 2003). In appendixIBlwe 
comment on testing the assumption of the statistical indepen- 
dence of a simulation patch. The second assumption expresses 
the notion of homogeneity. Again, in a real star, rotation or 
large-scale magnetic fields might lead to deviations from ho- 
mogeneity. The last assumption is a statement about ergodicity, 
which appears justified considering the chaotic nature of stellar 
granular flows (Steffen & Frevtaa 1995). All in all, we feel that 
the above assumptions are not too critical in view of the preci- 
sion we hope to achieve. Instead limitations of the simulations 
themselves - in particular the limited statistics due to the finite 
length of the calculated time series - control the accuracy of 
predictions at the moment. 

3. The radiative output of a simulation patch 

In the hydrodynamical simulations, the radiative energy ex- 
change is treated in great detail: the geometry of the radiat- 
ing flow is taken into account by a ray-tracing technique, and 
the wavelength dependence by a number of wavelength bands 
representing the effects of radiative heating or cooling in the 
continuum and spectral lines. While we concentrate in the fol- 
lowing on wavelength integrated ("white light") quantities, the 
formulas derived below also hold for variables which are mea- 
sured at a particular wavelength, or which are integrals over a 
wavelength interval. 

The basic radiative output that a simulation provides and 
that is relevant here, are surface averages of the emergent in- 
tensity I m for various inclination cosines fj. m , m e {1, . . . , M); 
i.e. the detailed spatial intensity information of a simulation 
is reduced to a level familiar from ID standard plane-parallel 
model atmospheres. However, the average intensities are com- 
puted for every time step, and are thus provided as a function 
of time I m (t). 

The discrete set of inclinations and azimuthal angles is 
rather small reflecting limits imposed by the available comput- 
ing resources. Typically, we worked with two or three angles in 
H~ and four in (^-direction. Table ^ summarizes the employed 
inclinations and associated weighting factors for evaluating the 
radiative flux. The integration scheme with, in total, three in- 
clinat ions M — 3 is a Gauss-Lob atto type integration scheme 
(e.g., lAbramowitz & Stegunll 19721) . while M = 2 is not. The 
M — 2 case is computationally quite economic for represent- 
ing the center-to-limb behavior in the flux calculation but has 
some disadvantages when it comes to calculating the photocen- 
tric variability, as we discuss in detail later. 

4. The statistics of the emergent flux of an 
individual simulation patch 

In the following, we derive expressions both for the tempo- 
ral average of the emergent flux of a simulation patch and for 
its standard deviation. Within the discretization of the radiative 
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Table 1. Inclination cosines //,„ and weighting factors w m for 
the flux integration schemes with a total number of inclinations 
of M = 2 and M = 3. Factor w m /if n is given due to its impor- 
tance for evaluating some disk-integrated quantities. 





M 


= 2 




M = 3 




m 


1 


2 


1 


2 


3 


Mm 


1.0000 


0.3333 


1.0000 


0.7651 


0.2852 


w m 


0.2500 


0.7500 


0.0667 


0.3785 


0.5549 




0.2500 


0.0833 


0.0667 


0.2216 


0.0451 



and 



<7j. = 07 V k,m 



(6) 



where cr s denotes the standard deviation of a variable g accord- 
ing to 



(7) 



Using the isotropy, we finally obtain for the temporal average 
of the emergent flux 



transfer equation that is used in the hydrodynamical simula- 
tions, the horizontal average of the emergent flux F (more pre- 
cisely of its vertical component) of a simulation patch is given 

by 



^ ^1 fin M N m 

= (hdClfil = I dfi I dipfil ^2n V V w km /j. m I k , 

J Jo Jo 



(1) 



where Wkm is the fraction of the solid angle of the hemisphere 
at inclination //,„ and at azimuthal angle <pi (see Fig. and 
h m the intensity in that direction. Here and in the following, 
the symbol denotes the approximate equality between con- 
tinuous expressions and discrete analogs, where the induced 
errors vanish in the limit of infinitely fine discretization. We 
denote physical approximations by Statistically, the con- 
vection pattern is horizontally isotropic. Therefore, all weights 
belonging to the same //-ring m (for the definition see Fig.Q are 
made equal since the discretization in the azimuthal direction 
should reflect the isotropy. The weights differ between differ- 
ent //-rings. Introducing the total weight w m of the //-ring m, the 
weights Wkm are normalized so that 



M N,„ 
m=\ k=l 



^ N m W km = 2^ Wm = 1 



(2) 



where N m = w m lwk m is the total number of solid angle ele- 
ments belonging to the //-ring m. For later reference, we give 
the discrete analog of the integral operator associated with the 
//-integration 



^1 M 

I dfxgfji) - V 

Jo ±T, 



W m g„ 



(3) 



where g is an arbitrary function of // and g m = g(ju m ). With Q 
we obtain for the temporal average of the radiative flux F of a 
patch 



M N,„ 

(F) =a 2n W kmMm (hm) , 

m=l k=\ 



(4) 



since w^m and //,„ are time-independent. Here and in the follow- 
ing, (g) denotes the temporal expectation value of a quantity g. 
Due to the horizontal isotropy, all elements of the solid angle 
belonging to the same //-ring share the same statistics 



WmjJ. m (I m ) . 



(8) 



Now, we want to obtain an expression for the expectation 
value of F 2 introducing variances and correlation coefficients 
among the intensities going into different directions of the solid 
angle. The variances and correlation coefficients could be eval- 
uated from the data provided by the hydrodynamical simula- 
tion. Of course, if one were interested in the flux variations of 
the simulation alone, one could resort to the flux directly eval- 
uated from Eq. Q and would not - rather indirectly - turn to 
the variances and correlation coefficients of the individual in- 
tensities. We derive the equations to compare them later with 
expressions for quantities integrated over the stellar disk. Using 
the basic relation Q we obtain 

M M N m N m i 

IF 2 } =i (27T) 2 ^ ^ 2j XI W k>n W k>m>HmlJ;n> ihmh'm') ■ (9) 
m=l m'=l k=l k' = l 

We introduce the linear correlation coefficient between I km and 
hm' C [hm, h'm'\ obeying the relation 

ihmh'm') - {hm){h'm'} + 0-i lm (Ti k , m ,C[hm,h'm>] 

= Qm )(Im') + C7„, C7 m , C[hmJk'm>] (10) 

with which we can rewrite expression as 

M M N m N m i 



m= 1 m'=l k=\ k' = \ 

+ O" /,„ 0" /„,/ C [hm, h'm'fj 

M M N„, N,„> 

= (F) 2 + (2n) 2 J]J]22 



m-l m' = l k=l k'=l 
WkmWk'm'UmUm' 0~ I,„ V /,„, C [hm, h'm' ] ■ 



(ID 



(hm) = (Im) V k, m 



(5) 



To progress on the analytical side, we consider two limiting and 
one intermediate case for the behavior among the individual in- 
tensities emerging from a simulation patch: 

i) The intensities of all elements of the solid angle are uncorre- 
lated, i.e., 

C [hm , h'm' ] = 5mm' Skk' Vm, Ttl' , k, k' . (12) 

ii) All intensities are perfectly positively correlated, i.e., 
C[/ ta ,W] = l Vm,m'XK. (13) 
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5.1. General statistical relations 




Fig. 1. Illustration of the weighting entering the calculation of 
the radiative flux. The plot shows the projection of a hemi- 
sphere, the concentric rings (in the text referred to as "yu-rings") 
and radial lines delineate the discrete elements of the solid an- 
gle (respective surface area elements for the integration over 
the stellar disk) entering the flux integration. Crosses mark the 
inclinations used for representing the angular dependence of 
the radiation field. The rectangles show how a square-shaped 
patch positioned on the surface at the various inclination an- 
gles appears foreshortened along the line of sight towards the 
limb. The diagram depicts the case M — 3 of Tabled 



iii) The intensities are correlated for all azimuthal angles be- 
longing to the same inclination fi, but are uncorrected among 
different inclinations, i.e., 



C [/fen, h'm'] = <W Vm, ml , k, k'. 



(14) 



Case ii) often provides a good description of the actual behavior 
found for the intensities of an individual patch. Case iii) is also 
a plausible behavior, since different inclinations correspond to 
different heights in the atmosphere; one would expect that first 
the correlation between different height levels is lost before the 
azimuthal correlation is reduced. In the three cases we obtain 
for the variance of the flux of an individual patch 



{2k) 2 



JiB=l "in W m^m cr /,„ 



(15) 



if C[hm, 

if C [I km , 
if C [I km , 

According to dl 11 1, o> is the sum of positive quantities so that 
as relation among the cases we find 



t 1k'm' ] 
t 1k'm' ] 
t 1k'm' ] 



(16) 



reflecting the higher degree of cancellation for increasing de- 
gree of de-correlation. 

5. Statistics of the disk-integrated, observable 
quantities 

We now want to derive expressions for the brightness and po- 
sitional fluctuations that a distant observer would record when 
observing the unresolved stellar disk. 



We are confronted with two problems: first, we have to intro- 
duce a correlation length on the surface of the sphere describ- 
ing the two-point spatial correlation of the temporal granular 
intensity signal. This we have done - as outlined before - by 
considering the stellar surface as composed of identical patches 
which are statistically independent from one another but share 
similar statistical properties. We replace the gradual behavior 
of the correlation function by a step-like behavior which sim- 
plifies the treatment of the problem tremendously. Second, we 
have to make best contact to the data that the simulations pro- 
vide. The central point is how we describe the center-to-limb 
variation. The basic idea is to transfer the discretization of the 
solid angle used in the hydrodynamical simulation to a corre- 
sponding partitioning of the stellar disk; i.e., we interpret the 
discretization of the solid angle from Table ^ as a partition- 
ing of the visible stellar surface, as depicted in Fig.^ This is 
of course analogous to the way the flux of a standard plane- 
parallel atmosphere model is translated to the disk integrated 
flux. Since we use a similar discretization of the stellar surface 
as in the case of integration over solid angle, the statistical for- 
mulae describing the flux of an individual simulation patch and 
disk-integrated emission will look almost identical. This is in- 
tentional. However, as will become apparent below, there are 
a number of subtleties which deserve some attention and may 
affect the final outcome substantially. 

Algebraically, we proceed by starting with the continuous 
description of a quantity, and discretize into patches to intro- 
duce the finite correlation length. We then return to the contin- 
uous description, since the analytical manipulations are most 
easily done in this representation. To finally establish con- 
tact with the simulation data, we discretize again applying the 
scheme used in the hydrodynamical model. 

We start by listing the discrete and continuous versions of 
the involved integral operators. Using the patches as discrete 
elements, the surface integral of a function g(fi = cos §, <p) over 
the stellar hemisphere can be written as 

j)dAgQ i ,#) = ^dfi^d?R 2 g(iJ,§)^Y^Ag k (17) 

where gk = g(jik,tfk), an d MA an d fk the coordinates of 
patch k. The star's radius is R, and N is the total number of 
patches on the curved surface of the visible stellar hemisphere; 
all patches have the same surface area A = R 2 Afi A<p, so that 



NA = 2nR 2 . 



(18) 



For formal reasons we relate the surface tiling to a spherical 
coordinate system, but in fact do not envision the stellar sur- 
face tiled by segments suggested by such a system. The tiling 
should be considered as locally close to isotropic — as would 
be the case for the tiling of a plane by squares. We repeat that 
an explicit construction of a tiling turns out to be unnecessary. 

We now formally define intensity-weighted integrals J g of 
an arbitrary function g as 



J g = (hdAfilg. 



(19) 
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This definition is introduced to streamline our nomenclature. 
We shall later set g — 1, and g — x where x is the Cartesian x- 
position of a patch on the stellar disk. J\ is related to the stellar 
flux, and J x to the position of the stellar photocenter. 

Similar to the angular integration we now make assump- 
tions about the statistics of the patches. To preserve axisymmet- 
ric conditions, the temporal expectation of the intensity should 
be a function of fi alone 



</>0u, ?) = </>(//)■ 



(20) 



Analogously, the standard deviation of the intensity should 
obey 

(Ti(jl,ip) = (TiQj.), (21) 



5.1 .1 . Temporal expectation value of J g 

In the following we consider a time-independent function 
g(fi, <p). With the previous assumptions about the //-dependence 
of the intensity and its variance, we obtain for the temporal ex- 
pectation value of J g 

XI r*7.7l r* 1 

dfi J o d V R 2 n </> g = 2nR 2 J o dfi fi(I)g, (22) 

where the last equality holds if g is also (^-independent. 



5.1 .2. Temporal variance of J g 

The temporal expectation of Jg depends on the properties of the 
spatial correlation of the intensity, similar to the flux of a simu- 
lation patch as expressed by Eq. We now introduce the dis- 
cretization into patches d!7l >. which allows us to express these 
properties easily. The accuracy of this rests on the assump- 
tion that the patch size is small enough in comparison to the 
global stellar scale that position-dependent quantities (like the 
time-averaged intensity and g) vary little across the patch, but 
large enough that spatial correlations among statistical quanti- 
ties across patches are negligible. We obtain 

N N 

(jf) * A 2 J^Y J f 1 Wk>(hh')gkgk> 

k=\ k'=\ 
2 N N 

= {j g ) + A 2 ^^fi k fi k ,cr Ik cr h ,g k g k ,C[I k ,I k >], (23) 

k=\ k' = l 

where the intensity I k is the surface average over patch k in 
the direction towards the observer. In the present context sum- 
mation index k does not refer to a particular coordinate direc- 
tion but merely numbers the patches located on the stellar sur- 
face. This reflects the fact that we discuss correlations among 
the intensities emitted by different patches and not within one 
patch as considered in Sect. @] The discretization into patches 
is not related to the partitioning of the stellar disk as depicted in 
Fig-E This is introduced at a later stage. We are mostly inter- 
ested in the case where the individual patches are statistically 
independent, i.e., C[I k ,I k >] = 8 kk >. However, to enable com- 
parison with the formulae derived for the flux of an individual 



patch, we also consider the case of perfect positive correlation, 
i.e., C [I k , I k >] = 1. For the case of statistical independence, we 
obtain 

cr) ^A 2 Yfi 2 k cr\g 2 k ^AR 2 dfi d V fi 2 cr 2 lg 2 . (24) 
k~\ Jo Jo 

Note, that relation i24i becomes asymptotically exact if the 
spatial correlation length of cr 2 goes to zero. For a ip- 
independent g, we get 



cr 



2 = 2nR 2 A f dfifi 2 cr 2 g 2 = (2n) 2 R A N- 1 
Jo 



I 



dfif/ajg 2 



M 



{2n) 2 R A N- 1 



Z2 2 2 
w m^ m cr Im g m 



(25) 



m= 1 



where in the last equality we introduced the segmentation l|3}. 

For the second case of a coherent behavior of the patches, 
again specializing for the case of a ^-independent g, we obtain 

,2 



cr, 



: \ 2 f m 

2 D 4 



Z 



W m /J. m O~ nl g„ 



(26) 



In summary, for g(p., <p) = g(fi) in the two cases we obtained 

f AT 1 fidfi^ojg 2 if C [/*,/*] = 8 kV 
o) =(2n) 2 R 4 \ /rl n2 (27) 

[^dutwjgj if C[4,4<] = 1, 

or the discrete analog using (|3j 



cr) ^(2n) 2 R 4 




N~ Z„=i y»ml4n 8m if C [h, h>l = 6 kk , 

I W m fX m CT[ u g m J if C [I k , It] = 1 . 



(28) 



At this point we would also like to emphasize a general 
feature of the formulae derived here and in the following: one 
might (erroneously) conclude from formula d2 8t > that cr 2 . de- 
pends on the particular, and to some extent arbitrary, geometry 
of the simulation patch since jV and <Ti m depend on the patch 
size. However, only their combination cr^J V^V enters the re- 
lations. This combination is invariant to changes in the patch 
geometry, as long as our basic assumption holds that a patch 
represents a statistically independent surface element. Hence, 
cr 2 . is in fact independent of the particular patch geometry. 



5.1 .3. Temporal power spectrum of J g 

To characterize the statistics of J g further by its temporal power 
spectrum, we only consider the case of statistically independent 
patches. In the following, y(ai) denotes the temporal Fourier 
component of a time-dependent variable y(t), and y* the conju- 
gate complex. Taking the Fourier transform of both sides of the 
basic definition dl9> . we get for an Fourier component of J g 



£ N 

J g = (pdAfilg - YjAvAgk 

J k=\ 



(29) 



where we assumed that g is time-independent. Accordingly, for 
a frequency component of the power spectrum we obtain 



N N 



(30) 



k=l k' = \ 
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This expresses the power of a frequency component for a par- 
ticular realization of a convection pattern and its evolution. We 
now have to average over all possible realizations to estimate 
the expectation value of the power. Denoting the ensemble av- 
erage as (g) ens , we obtain 



N N 



(31) 



fc=i k'=i 



All realizations share the property that the amplitude of a par- 
ticular Fourier component 4 is the same. Among the realiza- 
tions, the 4-components only differ by a uniformly distributed 
random phase factor exp(i<pw). The difference between two 
uniformly distributed random phases is again a uniformly dis- 
tributed random phase. It follows that 



(tJl) =\h\\ll\(e m,k ') =lJ* k 6kk' 

\ K I ens 1 1 1 K 1 \ / ens K 



(32) 



The last equality comes about because for k + k' the ensemble 
average {e tb '^ w ^ corresponds to an average over the unit cir- 
cle in the complex plane which is zero. For k — k' the phase dif- 
ference is not random but always zero, and the ensemble aver- 
age is one. From now on we drop the subscript "ens" indicating 
the ensemble average, but keep in mind that with respect to av- 
erage power, angular brackets denote ensemble averages. One 
should further keep in mind that the hydrodynamical models 
provide only an estimate of the average power of the intensity 
J*/?, because only a finite time interval is simulated. 

With the formula for the ensemble average of the power of 
the intensity we obtain for the power of J g 



M) - J]Ay k (ij:) ; 

k=l 

- AR 2 J "d<pfi 2 (ll*)g 2 . (33) 

Similar to Eq. $241 the above relation becomes asymptotically 
exact, if the spatial correlation length of the intensity goes to 
zero. Assuming a ^-independent g, we finally obtain 

(j g J' g ) = Qxf&N- 1 j ' dpi? (It)? 

M 

- dnf^N- 1 J] *W4 (U* m ) gl- (34) 

m=l 

Note, that the above relation for the expectation value of a fre- 
quency component of the power spectrum of J g is consistent 
with the relation for its variance d28i : the frequency integral of 
the power (excluding zero frequency) equals the variance of the 
signal as expressed by Parseval's theorem. 

5.2. Statistics of the observable flux 

We now apply the previously derived formulae to evaluate sta- 
tistical properties of the observable flux. The flux / observed at 
a distance D from the star is related to the moment J\ according 
to J\ = D 2 f. Due to this close connection we shall sometimes 
loosely refer to J\ as observable flux. 



5.2.1 . The expectation value of the observable flux 

From Eq. \22\ setting g — 1 and the basic definition of the flux 
of an individual patch ([0, we obtain for a (^-independent mean 
intensity (I m ) 



(J 1 )=2ttR 2 f dnii(I)=R 2 (F}.. 
Jo 



or 



(35) 



(36) 



The observed temporally averaged flux - perhaps no sur- 
prise - is the average flux of an individual patch reduced by 
the distance-related factor R 2 /D 2 . At this point one might be 
tempted to ignore the fact that Eq. d36i expresses a relation 
between the expectation values of / and F, not the fluxes 
themselves, and might calculate a power spectrum of the disk- 
integrated flux as a power spectrum of the flux of the simulation 
patch, perhaps scaled by the factor A^ 1 . We discuss the quanti- 
tative error of this approximate procedure in Sect. [6] 

5.2.2. The variance of the observable flux 

To compute the variance of the observable brightness fluctua- 
tions, we need to obtain the second moment of Ji . As in the 
case of an individual simulation patch, we consider the two 
limiting cases of statistical independence and perfect positive 
correlation. Applying Eq. d28l > results in 



~ {2n) n 1 (vt* \ 2 -f ru i i 1 (37) 

We find - perhaps again no surprise - that in the case of sta- 
tistical independence the variance of the observable flux scales 
inversely proportional to the number of patches on the visible 
stellar surface. Relations J37I are close to the results for an in- 
dividual patch |15i . We can make it even closer by introducing 
the number of patches N m located within /i-ring m of the stellar 
disk 



N m 



w m N, 



(38) 



since all patches on the stellar surface have the same area. 

Since the weights in the solid angle integration are not 
proportional to the number of solid angle elements, Eq. J 1 51 
cannot, on the other hand, be cast in the form of Eq. i37\ . 
Introducing (1381 makes Eq. d37l > read 



<^-<W(§)' 



Zli N m x w 2 m ii 2 m cr 2 m if C [4, /*] = So. 

[Lm=i W m ll m CTi m j if C [4, 4<] = 1. 



(39) 



Only in the (for convective flows) unrealistic situation of a per- 
fect correlation among the intensities of a patch and also among 
the intensities across the stellar disk, one can directly relate 
disk-integrated brightness fluctuations to the fluctuations of the 
flux of the patch as crj, = R 2 cr F which is analogous to the re- 
lation for the flux {J\) = R 2 (F). In all other cases one has to 
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resort to the intensities themselves. Even if one chose to make 
all elements of the angular integration the same size, and ob- 
tained a direct correspondence of the formulae this way, one 
still could not use the flux of the patch to derive the brightness 
fluctuation of the star as a whole. The reason is found in the 
different degrees of correlation among the intensities of a patch 
and across the disk. Within a patch the intensities going in vari- 
ous directions usually show a substantial degree of correlation, 
which is lacking for the intensities emerging from different lo- 
cations on the stellar disk. We discuss this aspect quantitatively 
in Sect.|6] 

5.2.3. Temporal power spectrum of the observable flux 

To further characterize the statistics of the observable flux, we 
use Eq. J34i to write down its temporal power spectrum for the 
case of statistically independent patches 



(40) 



For the power spectrum of the observable, relative flux varia- 
tions we obtain 



(//*) _ (fity . 111^/4(4/;) 



(fr {JxY 



(Zm=l W mf*m ihn)) 



(41) 



We remind the reader that the angular brackets around fre- 
quency components of power spectra denote ensemble aver- 
ages. 

5.3. Statistics of the photocentric position 

We consider the position of the stellar photocenter in Cartesian 
coordinates (x,y). The origin of the Cartesian coordinate sys- 
tem is placed in the geometrical center of the stellar disk. We 
again assume statistically isotropic conditions with respect to 
disk center, x- and y-direction are equivalent, and we only have 
to regard - without loss of generality - one Cartesian coordi- 
nate, say x. The position of the photocenter x p h(?) is given as the 
intensity-weighted mean of the x-positions of all light emitting 
patches tiling the visible stellar surface according to 



1 +0 



<•/!> 



(42) 



Since the relative fluctuations of the observable flux are small, 
we neglect the higher order terms in the fluctuations of J\ and 
work with the second equality in Eq. \A2\ . This leaves us with a 
much more convenient linear relationship between x p h and the 
intensities. For obtaining the expectation value and variance of 
x p h we now have to evaluate the respective quantities of J x . 
We do this only for the case of uncorrelated patches, since it is 
clear that one would not obtain any motion of photocenter for 
the case of a perfectly correlated action of all patches. 

To apply the previously derived formulae, we set g = x, in 
spherical coordinates 



X = Ryf 



1 - fl 2 cos tp. 



(43) 



5.3.1 . Expectation value of the photocentric position 
From Eq. ( 1221 one immediately gets 

,-) = J dfi j^dtp R 3 ft (I) yjl -ii 1 cos (p = 0, (44) 



</ 
and 

(jCph) = 0, 



(45) 



since the ^-integral vanishes. The expectation value of the pho- 
tocenter's position coincides with the geometric center of the 
star. This is of course a consequence of the fact that we - in a 
statistical sense - preserved axisymmetry in the formulation of 
the problem. 

5.3.2. The variance of the photocentric position 

To obtain the variance of the photocentric position, we have to 
evaluate From relation (1241 . valid for ^-dependent quan- 
tities, we get 



cr] = AR* 



Jr>l pin 
dji I di 
o Jo 



dip f/oj ( 1 - y 1 ) cos z (p 



) h 



where we used the relation (see Eq. 1251 ) 



(InfR^N^ 



f 

Jo 



j 2 2 
djjfl (Tj. 



(46) 



(47) 



Equation ( I46> describes a relation between J x and J\ . We shall 
see later that the connection is even tighter than it might at first 
appear from Eq. (I46> . since it turns out that the term in angu- 
lar brackets is only weakly variable with the degree of limb- 
darkening of cr]. Now we can write the variance of x p h: 



2 2 

2 <r] 



(JiY 



4n 



f W'>) 



^dnn A cr] 



2^2 



(48) 



Remarkably, the variance of the photocentric position is 
independent of stellar radius (or equivalently total number of 
patches), as long as we assume that the patch size A and its ra- 
diation properties are fixed. Astrophysically, this means that a 
star's photocentric variability 1 depends only on its atmospheric 
parameters (T^jf, logg, and chemical composition). However, 
when changing a star's radius at fixed surface conditions by al- 
tering its mass, the relative importance of sphericity is changed. 
This is a higher order effect but might bring a radius depen- 
dence back in for objects with extended atmospheres. The ra- 
dius independence is intrinsic to the particular 2D axisymmet- 
ric geometry of the problem. It can be shown that the same 
does not hold, e.g., for a linear chain of patches. From a prac- 
tical point of view it is convenient that we do not have have to 



1 Note that we are discussing the absolute spatial variability of the 
photocenter, hence the stellar distance does not enter at this point. 
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Fig. 2. Comparison of the integral LdfiL 2 fi 2 , which enters the 
evaluation of cr^ (see Eq. 1471 ). For a typical limb-darkening 
around a = 0.6, we obtain a difference of 13 % to the exact 
result in the case M = 2. In the case M — 3, differences are 
negligible. 
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Fig. 3. Comparison of J^dpL 2 p 4 / J dfxL 2 fi 2 which enters the 
computation of the displacement of the photocenter (see 
Eq. 1481 ). For a typical limb-darkening around a = 0.6, we 
obtain a difference of 30 % to the exact result in the M = 2 
case. In the M = 3 case, differences are negligible. 



specify a precise stellar radius when evaluating the variance of 
the displacement of the photocenter. 

We have not written down the discrete analogs of Eqs. d46b 
and J48t > yet, because a numerical problem exists when evalu- 
ating integrals involving higher powers of p. The integrals be- 
come inaccurate; this is particularly the case for our integration 
scheme with M — 2, which is not a Gauss-Lobatto type integra- 
tion. We now seek a way to mitigate the problem. To illustrate 
the situation, we assume that the dependence of the expecta- 
tion value of the intensity on p follows a linear limb-darkening 
law L of the form 



COG") 



L(p) — 1 — a + a/j 



(49) 



</>0"=D 

where a is the limb-darkening coefficient. The numerical re- 
sults of the hydrodynamical simulation show that the standard 
deviation of the intensity is often roughly proportional to the 
intensity. Hence, it is reasonable to assume a similar behavior 
for the standard deviation 

07 (p) 



(50) 



o-iin = 1) 

Using the linear limb-darkening laws ( I49> and (I50> . the inte- 
grals in Eq. (I48> can be evaluated analytically (see appendixlAl. 
as well as numerically. The results can be compared to estimate 
the error introduced by the various integration schemes. 

Figure |2 shows the comparison between the exact and nu- 
merical results for the integral J () dpL 2 p 2 , which enters the cal- 
culation of a 2 (Eq. 1471 ). For a typical limb-darkening at op- 
tical wavelength we find a sizeable but tolerable integration er- 
ror in the case M — 2. Case M-3 gives almost exact results. 
Consequently, we do not expect that the limited accuracy of our 
integration schemes will affect the evaluation of the photomet- 
ric variability substantially. 

As shown in Fig. the situation deteriorates when con- 
sidering the ratio Jj^LV/ J Q dpL 2 p 2 , which enters the eval- 



uation of the displacement of the photocenter (Eq. 1481 ). 
However, the figure also shows that the exactly evaluated over- 
all variation of the ratio of the two integrals is not large. In view 
of the large error associated with the numerical integration in 
the case M = 2, this leads to the idea of replacing the ratio 
of the integrals by a constant, representative value of the exact 
integral. We choose for the ratio 



2 
3' 



(51) 



which results from a limb-darkening coefficient of about a = 
0.6. Apart from a = 0.6 being a reasonable value for the limb- 
darkening coefficient in the optical, a ratio of | is chosen as 
a reminder that the approximation i5 11 is of course closely 
related to the well-known Eddington approximation, where a 
value of | is assumed for the ratio between the second and ze- 
roth angular moment of the intensity. For this reason we shall 
refer to approximation < I5 1 1 as the Eddington-like approxima- 
tion. 

We apply the Eddington-like approximation to relations for 
the displacement of the photocenter J46> and J48ft . In doing so 
we hope to mitigate the effects of decreased numerical accu- 
racy when evaluating the ratio J Q dpL 2 p 4 / J Q dpL 2 p 2 . For the 
standard deviation cr Jx , we obtain 



07, 



R 



07, 



(52) 



or for the standard deviation of the photocenter cr x h itself, 



(53) 



^ * y/6Vi)~ V6</>' 

Equation (I53> relates the observable flux variations to the dis- 
placement of the photocenter. Remarkably, no information in- 
trinsic to granulation enters this relation. Hence, the relation 
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could be more generally valid t han the particula r derivation 
might suggest. And indeed, e.g., Lindegre nl dl977t) has derived 
a similar connection in the context of the astrometry of objects 
with an apparently time-variable surface pattern due to rotation. 
Equation ( I53l l is a relation among quantities which should be 
directly accessible to observation. This might be considered an 
advantage. However, concerning convection as such, the diag- 
nostic power of Eq. (1531 is limited since the observed variations 
in brightness and position cannot be immediately connected to 
convection. Other processes can produce a similar signal. 

Since the temporal power spectrum of the brightness fluc- 
tuations is closely linked to its variance, Eq. d53l also states 
that the power spectrum of the photocentric displacement is - 
up to the frequency independent factor R 2 /6 - identical to the 
power spectrum of the relative brightness fluctuations. This is 
a very convenient property. The similarity between the power 
spectra of brightness and positional fluctuations does not imply 
a temporal correlation among these variables. In fact, one can 
immediately show that C [J\ , J x ] = or, correspondingly, that 
the observable flux and the displacement of the photocenter are 
uncorrelated: 

C[/,* ph ] = 0. (54) 

lHatzesl J2002) discusses simultaneous measurements of corre- 
lated observables to improve signal-to-noise ratios in exoplanet 
searches. The observables discussed here do not offer such a 
possibility. 

Using the Eddington-like approximation in Eq. d48i . we al- 
ternatively obtain 



d/2 fx cr 



fT2n jT^ M (/> 



_1_ 

(E,1f = i"w>i) 2 

■ 2m=l W m fl m (I,„) 



(55) 



where we have introduced the linear patch size I = Va. A gain, 
formula d55i is manifestly radius-independent. It also gives the 
discrete expression for evaluating the standard deviation of the 
photocentric displacement. 

6. Estimating cr f from cr F ? 

In this section we briefly digress to discuss the relation between 
the flux variations of a simulation patch as given by Eq. dl5> . 
and the disk-integrated flux variations given by Eq. (I37> more 
quantitatively. The reason is that one might take the standard 
deviation of the flux of a simulation patch scaled by 1 / V^V 
as a proxy for the standard deviation of the observable flux. 
One might proceed in the same way when calculating power 
spectra of the disk-integrated flux variations. For deriving an 
estimate of the introduced errors, we once again assume the 
linear limb-darkening laws ( 1491 and ( 15 0> . We evaluate the ratio 
of the relative standard deviations scaled by the total number 
of patches VAf(cry//)/(cr F /,F). We consider the cases ii) and 
iii) (see Eq. 1151 ) of the angular correlation of the intensities 
of a patch, and compare those with the (realistic) case of a sta- 
tistically independent behavior of the patches across the stellar 
surface. We restrict ourselves to cases ii) and iii), because they 
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Fig. 4. The disk-to-patch ratio of relative flux variations as a 
function of the limb-darkening coefficient for two different as- 
sumptions about the angular correlation of the intensity within 
a patch. 



encompass the actual range found for the angular intensity cor- 
relations within a patch. Figure 0] depicts the result. The exact 
result, i.e., analytically integrated angular integrals, is left out 
in the case of partial correlation, since no continuous analog of 
the discrete formula d 1 51 (case iii) can be easily written down. 
It is of minor importance since what is relevant here are not the 
deviations from the exact result, but the absolute value of the 
disk-to-patch ratios. 

Figure |4] shows that one tends to underestimate the ampli- 
tude of the observable brightness fluctuations, if one computes 
the fluctuations directly from the flux of a simulation patch and 
simply scales it by 1 / y/N. Depending on the angular corre- 
lation of the intensities of the simulation patch, the error can 
become substantial. In the (not shown) case of completely un- 
correlated intensities of a simulation patch, the error becomes 
even larger (a ratio of about 2 for M — 2, about 3.3 for M = 3, 
both for four azimuthal angles) than in the cases depicted in 
Fig-E] The two reasons for the discrepancies have already been 
discussed in connection with Eqs. d37t and fl39l >. 

7. An analytical model for the case of a linear 
limb-darkening law 

It is instructive to evaluate the formulae derived in Sects. 15.21 
and 15.31 for the linear limb-darkening laws for / and 07 
(Eqs. d49t and 1501 ). From Eqs. i22\ . d28i . and auxiliary formu- 
lae from appendixEl we directly obtain for the relative fluctua- 
tions in the observable flux for the case of uncorrelated patches 



cr f 1 cr, 

— = N~~- — 
if) (I) 



& L2 4 _x err VFI^i 



- err 



(56) 



Here (/) and cr; should be understood as the values at /1 = 1 . 
The a-dependent numerical factor in Eq. d56t has already been 
depicted in Fig.0] the solid curve depicting the exact integral 
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in the case of complete correlation. Dependence on the limb- 
darkening coefficient a is only mild, but it shows that the bright- 
ness variations grow for a more pronounced limb-darkening. 

Similarly, we obtain from Eq. J48I for the standard devia- 
tion of the photocentric displacement 



1 

(jjjWaV)) 5 



' A|, h 



foul? 

\l_ _ a_ ~ £_ 

I CT[ y 15 10 + 42 



(57) 



The fl-dependent factor in Eq. (I57> is only very weakly vari- 
able, as it varies (non-monotonically) between 0.201 and 0.206 
with the limb-darkening coefficient in the interval a e [0, 1] so 
that to good approximation 



0V„ » 021%. 
ph (I) 



(58) 



Formula d58l > clearly emphasizes that the photocentric displace- 
ment does not depend on the stellar radius but only on local, 
granulation-related properties. 

The chosen box-geometry of the hydrodynamical simula- 
tions sets the linear patch size, and the patch-averaged tem- 
poral intensity fluctuations scale with the surface area of the 
box. The natural statistically independent elements on the stel- 
lar surface are not patches but the granules, so that in phys- 
ical terms one would like to connect / to the linear size of a 
granule and 07/ (/) to the time-wise fluctuations of its emit- 
ted light. In the spirit of our basic assumption of ergodicity, we 
might further replace the time-wise fluctuations by equivalent 
space-wise fluctuations. We preliminarily find an approximate, 
empirical relation 



/£i*04/ ^2 

' ~ v.t (gran , 



which, combined with Eq. J58l >. gives 



O-.v 



0.08 L 



SI n 



(59) 



(60) 



Here Z glan is the linear granular size - defined via the maxi- 
mum in the spatial power spectrum of the intensity pattern - I 
the spatial average of the intensity, and d>/ lms the spatial root- 
mean-squares intensity contrast - (conventionally) defined at 
ideally infinite resolution, here the resolution of the compu- 
tational grid. Equation ( I59t approximately relates the patch- 
connected quantities in Eq. (I58l > to physical properties of in- 
dividual granules and the granulation pattern. Perhaps surpris- 
ing at first glance, the numerical factor in Eq. d59l is not one 
but smaller than one, despite the fact that the product Icrj / (/) 
should be scale-invariant, as long as no spatial correlations 
are present on the considered scales. The intensity contrast in 
Eq. ( 1591 refers to scales smaller than the granular scale where 
this assumption is violated. On these scales the intensities are 
spatially correlated leading to enhanced fluctuations relative to 



larger scales, which needs to be compensated for by a numer- 
ical factor smaller than one. Moreover, we also expect devi- 
ations from one due to the differences in geometry between 
square-shaped boxes and roundish granules. 

As a rough estimate for the Sun from Eq. i60\ with Z gran « 
1000km and SI lms /I ~ 0.18, we obtain a granulation in- 
duced photocentric variability of cr Aph ~ 14 km, * 10~ 7 AU, or 
« 0.1 yuas / D [pc]. This already indicates that granulation will 
hardly affect the astrometry of dwarf stars at the precision level 
that is expected to be reached by the Gaia mission. The size 
of granul es grows roughly in versely proportional to the surface 
gravity ( Frevta g et al. I ll997t) of a star, so that in giants and su- 
pergiants we expect much greater effects. 

8. Conclusions and outlook 

We derived formulae for evaluating the disk-integrated fluctu- 
ations of a star's brightness and the position of its photocenter 
from local-box hydrodynamical atmosphere models. Besides 
the merely technical results, we showed that there is a close 
connection between stellar photometric and photocentric vari- 
ability. While derived here for granular surface patterns, the re- 
lation should be similar for any statistically homogeneous time- 
variable surface pattern. Activity related brightness patterns are 
an interesting case. Unfortunately, active regions are not ex- 
pected to be homogeneously distributed over the stellar sur- 
face but to follow belt-like zones. We nevertheless expect that a 
semi-quantitative connection will still remain between activity- 
induced photometric and photocentric variability so that order- 
of-magnitude estimates of the photocentric variability are pos- 
sible from measured brightness fluctuations. This may help in 
estimating limitations to the achievable astrometric accuracy 
for missions like Gaia. As we indicated in the previous section, 
effects due to granulation themself are not likely to be impor- 
tant with the exception of extreme giants. This will be quan- 
tified in more detail in a subsequent pape r. A preview of th e 
simulation results can be found in Svensson & Ludwig (2005). 
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Appendix A: Angular integrals of L n /j. m 

For a linear limb-darkening law of the form 

L(jj) = 1 - a + an, (A.l) 
we find the following integrals over 

(A.2) 
(A.3) 
(A.4) 
(A.5) 
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Appendix B: Statistical independence of a patch 
and the spatial autocorrelation of the 
granulation pattern 

One of our basic assumption is that the radiative output of a simulated 
patch can be regarded as statistically independent of all other patches. 
We might check this assumption by testing whether the assumption is 
already fulfilled within a simulation patch, i.e., whether parts of the 
patch are already distant enough to radiate independently. For this we 
have to evaluate the time-wise correlation coefficient of the emitted 
intensity at two locations (x, y) and (x',y r ) given by 

C[I(x,y),I(x',y')] 

(I(x,y)I(x',y'))-(I(x,y))(I(x',y')} 

= = (g(x,y)g(x ,y )> , (B.l) 

0"/(*,v)0"/(A-'y) 

where we introduced the normalized intensity 
/-</} 

g = — . (B.2) 

07 

Due to the horizontal statistical homogeneity of the granulation pat- 
tern, the correlation is expected not to depend on the absolute position 
but only on the distance between two points. To get an overall measure 
of the statistical dependence C among the various locations, one can 
spatially average over all pairs of points located a distance (Ax, Ay) 
apart: 

C(Ax,Ay) 

= - i -If dx f dyg(x,y)g(x + Ax,y + Ay)).(B.3) 

O2 - xi)(y 2 -yi) \J Xi J V1 / 

C is given by the time-wise average of the autocorrelation function 
of the normalized granular intensity divided by the considered surface 
area. As stated before, C depends on distance alone so that one can 
further average the above expressions over circles Ax 2 + Ay 2 = const., 
and obtain a one-dimensional function describing the decrease of cor- 
relation with increasing distance between two locations. In order to 
have our basic assumption fulfilled, the spatial scale over which the 
correlation drops should be significantly smaller than the linear patch 
size. 

Here we have regarded the intensity, not the flux. However, since 
all inclinations follow roughly the same behavior, the correlation 
length of the intensities of all inclinations should be roughly the same. 
This implies a corresponding correlation length for the radiative flux. 
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